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Abstract 

We examine the long-wavelength current response in anisotropic superconductors and show how 
the field-dependence of the Meissner penetration length can be used to detect the structure of 
the order parameter. Nodes in the excitation gap lead to a nonlinear current-velocity constitutive 
equation at low temperatures which is distinct for each symmetry class of the order parameter. 
The effective Meissner penetration length is linear in H and exhibits a characteristic anisotropy for 
' fields in the a6-plane that is determined by the positions of the nodes in momentum space. The 

nonlinear current- velocity relation also leads to an intrinsic magnetic torque for in-plane fields that 
are not parallel to a nodal or antinodal direction. The torque scales as H 3 for T —* and has a 
characteristic angular dependence. We analyze the effects of thermal excitations, impurity scattering 
,.0 , and geometry on the current response of a d x 2_ y 2 superconductor, and discuss our results in light 

D ■ of recent measurements of the low-temperature penetration length and in-plane magnetization of 

single- crystals of YBa^CusOis and LuBaiCuzOi-s. 
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I. INTRODUCTION 



Recent measurements of the Meissner penetration depth [jl[ and Josephson interference effects in Y Ba2Cu^OT-s 
' H have been interpreted in support of a spin-singlet order parameter belonging to the one-dimensional, d x 2_ y 2 
ON ■ representation, A(py) = Ao — K,), which breaks reflection symmetry in the basal plane. Such a pairing state has 
been proposed by several authors ||-|| based on arguments that the CuO materials are Fermi liquids close to an SDW 
| instability. 

If the cuprates have an order parameter that is unconventional, i.e. one that breaks additional symmetries of 
the normal state besides gauge symmetry, then the superconducting state is expected to exhibit a number of novel 
properties, including (i) gapless excitations below T c , (ii) anomalous Josephson effects, (iii) exotic vortex structures 
and associated excitations, (iv) new collective modes, (v) sensitivity of superconducting coherence effects to defect 
scattering and (vi) multiple superconducting phases. j^|?J Many of these signatures of unconventional pairing have been 
observed in superfluid 3 He, and in heavy fermion superconductors, notably UPt3. Q The case for an unconventional 
. i-H , order parameter in the cuprates, and particularly a d x 2_ y 2 state, is not settled; there are conflicting interpretations of 
closely related experiments, variation in results that are presumably related to material quality or preparation, 

j pd| , ^2| and experimental results that are not easily accounted for within the d x 2_ y 2 model. [^3 14 



In this paper we examine the long-wavelength current response in superconductors with an unconventional order 
parameter, and show how the field- dependence of the Meissner penetration length can be used to detect the structure 
of the order parameter. This report extends our earlier work on nonlinear supercurrents, f p^|p^ ] and provides the 
relevant analysis that could not be included in our short reports. Specifically, we show (i) how the nodes in the 
excitation gap, whose multiplicity and position in momentum space depend on the symmetry class of the order 
parameter, lead to a nonlinear current-velocity constitutive equation at low temperatures (T <C T c ) which is unique 
and qualitatively distinct for each symmetry class. The effective Meissner penetration length is linear in H and 
exhibits a characteristic anisotropy for fields in the afe-plane. (ii) This anisotropy is determined by the positions of 
the nodes in momentum space. For example, in the case of a d x 2_ y 2 state in a tetragonal material the anisotropy is 
precisely l/y/2, independent of the detailed shape of the Fermi surface or the gap. (iii) The nonlinear current- velocity 
relation leads to an intrinsic magnetic torque for in-plane fields that are not parallel to a nodal or antinodal direction. 
The torque scales as H 3 for T <^T C and has a characteristic angular variation with a period of 7r/2 (for tetragonal 
symmetry). The magnitude and angular dependence of this torque are calculated for thick superconducting films or 
slabs, (iv) We discuss the effects of thermal excitations, impurity scattering and geometry for observing these features 
in a d x 2_ y 2 superconductor. Recent measurements of the low-temperature, zero-field penetration length [jj] are used 
to determine the relevant material parameters for YBazCusOi-s, which are then used to estimate the magnitudes of 
the field-dependence of the penetration depth and the torque anisotropy at low temperatures. 
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Our starting point is Fermi-liquid theory applied to anisotropic superconductors; section II includes the relevant 
theoretical framework needed to calculate the current response in unconventional superconductors. We derive formulas 
relating the equilibrium supercurrent to the magnetic field and discuss the linear response limit in section III. The 
nonlinear current-velocity constitutive equation is examined in section IV. A clean superconductor with a line of nodes 
in the gap has an anomalous contribution to the current which is a nonanalytic function of the condensate velocity 
v s , at T = 0. The relation of the anomalous current to the quasiparticle spectrum is discussed, and the contribution 
of this current to the Meissner penetration depth is obtained from solutions to the nonlinear London equation. 
The effects of impurity scattering and thermally excited quasiparticles on the anisotropy and field-dependence of the 
supercurrent are examined in detail; the signatures of the anomalous current survive thermal excitations and impurity 
scattering at sufficiently low temperatures and weak (or dilute) impurity scattering. We discuss our results in light of 
recent experiments on the low-temperature penetration depth Q in single crystals of YBa2Cu^OQ ^. An important 
conclusion is that if the linear temperature dependence of the penetration depth reported for YBa2Cu30§,$b is due 
to the nodes of a d x i_ v i order parameter, then the nonlinear Meissner effect, including the intrinsic anisotropy and 
field-dependence, should be observable for T < IK with a change in A a b of approximately 30 A over the field range 
< H < H c i ~ 200 G. In section V we discuss the nonlinear current, and associated in-plane magnetic torque, 
that develops for surface fields that are not aligned along a nodal or antinodal direction. The torque anisotropy 
(or transverse magnetization) is obtained from solutions to the nonlinear London equation at low temperatures. 
We also comment on a recent experimental report of a measurement of the in-plane magnetization fl7|| of a single 
crystal of LuBa2Cu^,Oj^s- In the rest of the introduction we briefly discuss the symmetry classes and unconventional 
order parameters for superconductors with tetragonal symmetry appropriate to the CuO superconductors (see Refs.( 
p8|,[7pS|]) for detailed discussions.) 

Symmetries of the pairing state 

BCS superconductivity is based on a macroscopically occupied equal-time pairing amplitude f a p{pj) ~ 
(a,f} fCl a-fiff)}, for quasiparticle pairs near the Fermi surface with zero total momentum and spin projections a and /3. 
Fermi statistics requires that the order parameter obey the anti-symmetry condition, f a p(pf) = —f/3 a (~Pf)j while 
inversion symmetry (if present) implies that the pairing amplitude decomposes into even-parity (spin-singlet) and 
odd-parity (spin-triplet) sectors. Furthermore, the pairing interaction separates into a sum over invariant bilinear 
products of basis functions for each irreducible representation of the point group. The resulting ground-state order 
parameter, barring the exceptional case of near degeneracy in two different channels, belongs to a single irreducible 
representation. For tetragonal symmetry there are four one-dimensional (ID) representations and one two-dimensional 
(2D) representation, and each of them occurs in both even- and odd-parity representations.^] The residual symmetry 
of the order parameter is just that of the basis functions for the ID representations, but for the 2D representation 
there are three possible ground states with different residual symmetry groups. There is no evidence that we are 
aware of to support a spin-triplet order parameter in the CuO superconductors; in fact the temperature dependence 
of the Knight shift in the cuprates ff(],|2l[] is argued to strongly favor a spin-singlet order parameter. |^2| Thus, we 
limit the discussion to even-parity, spin-singlet states; however, most of the analysis and many of the main results for 
the current response are also valid for odd-parity states. 



TABLE I. Even Parity Basis functions and Symmetry Classes for D/th 
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Table I summarizes the symmetry classes of the order parameter for spin-singlet pairing. All of the ID repre- 
sentations have residual symmetry groups which include four-fold rotations combined with appropriate elements of 
the gauge groups. The states E g (l,0) and E g (l, 1) have a residual symmetry group that allows only two-fold rota- 
tions. The resulting supercurrent, or superfluid density tensor, for such states is in general strongly anisotropic in the 
basal plane. The 2D order parameter, E g (l,i), preserves the four-fold rotational symmetry, but breaks time-reversal 
symmetry. 

Although the B\ g (d x 2_ y 2) and B2 g {d xy ) order parameters break the C4 rotational symmetry of the CuO planes, a 
combined C4 rotation and gauge transformation by e m is a symmetry. Since many properties of the superconducting 
state depend only on Fermi-surface averages of |A(p)-)| 2 , the broken rotational symmetry is not easy to observe. In 
particular, the London penetration depth tensor is cylindrically symmetric for any of the ID pairing states listed in 
Table I. Furthermore, all of the unconventional gaps in Table I yield a linear temperature dependence at T -C T c for 
the zero- field penetration depth (in the clean limit). 

A distinguishing feature of each phase, which is a consequence of their particular broken symmetries, is that the 
nodes of each gap are located in different positions in p-space. A point that we make below is that the field-dependence 
of the supercurrent may be used to locate the positions of the nodal lines (or points) of an unconventional gap in 
momentum space. This gap spectroscopy is possible at low temperatures, T <C T c , and is based on features which are 
intrinsic to nearly all unconventional BCS states in tetragonal or orthorhombic structures. 



II. FERMI-LIQUID THEORY OF SUPERCONDUCTIVITY 

Our starting point for calculations of the current response is the Fermi-liquid theory of superconductivity. This 
theory is general enough to include real materials effects of Fermi-surface anisotropy, impurity scattering and inelastic 
scattering from phonons and quasiparticles, in addition to unconventional pairing. A basic feature of the Fermi-liquid 
theory of superconductivity (c.f. Refs.( p3|-p5|) for a more detailed discussion of the formulation of Fermi- liquid 
theory.) is that for low excitation energies (hu>, ksT 1 hqvf, A) <C Ef, the wave nature of the quasiparticle excitations 
is unimportant and can be eliminated by integrating the full Matsubara Green's function over the quasiparticle 
momentum (or kinetic energy) in the low-energy band around the Fermi surface, 

g a p(Pf,R; en) = - r <% f dr e^ T ' h [ dfe~^ h < T T ^ a {R + r/2, r) ^(R - f/2, 0) > , (1) 

J-uj a Jo J 

where £p = Vf{pf)(\p\ — \pf\) is the normal-state quasiparticle excitation energy for momentum p nearest to the 
position pf on the Fermi surface and Vf(pf) is the quasiparticle velocity at the point pf. The resulting quasiclassical 
propogator is a function of the momentum direction pf on the Fermi surface, the center of mass coordinate R and the 
Matsubara energy e„ = (2n + 1)ttT. The pairing correlations are described by the ^-integrated anomalous Green's 
functions, 



f a p(p f ,R; e n ) = - I dtp f dr e^' h I df e -^ h < T T i> a {R + r/2, r) ^{R - f/2, 0) > . (2) 







The low-energy quasiparticle spectrum, combined with charge conservation and gauge invariance, allows one to formu- 
late observables in terms of the quasiclassical Green's function and material parameters defined on the Fermi surface. 
For example, the equilibrium current is given by 

l(R) = -eN f f dp f v f {p s )TY, \ Tr {f 3 g(p f , R; e n )} , (3) 

where Nf is the single-spin density of states at the Fermi level and the integration is over the Fermi surface with a 
weight factor of the angle-resolved density of states normalized to unity. We have introduced the 4x4 quasiclassical 
functions in 'spin x particle-hole' space; a convenient representation for the particle-hole and spin structure of the 
propagator is, 

~ = ( 9iPf,R\tn) + g(pf,R\_e n ) ■ 3 f(l>f, Rj e n ) io v + f(pf, Rj e„) • ido v \ ^ 

V R; e«)* ivy - f(-I>f,R; e„)* ■ ia y a g(-p f ,R; -e n ) - g(-pf,R; -e n ) ■ a y aa y J 

This matrix structure represents the remaining quantum mechanical degrees of freedom; the coherence of particle 
and hole states is contained in the off-diagonal elements in eq. ([!]). The diagonal components are separated into spin- 
scalar, g, and spin-vector, g, components. The scalar component determines the current response, while the vector 
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components determine the spin-paramagnetic response. The off-diagonal propagator separates into spin-singlet, /, 
and spin-triplet, /, pairing amplitudes, which are coupled to the diagonal propagators through the quasiclassical 
transport equation, 



Q[g , o-] = \ie n T 3 - a(pf, R; e n ) , g(p f , R; e n ) + iv f ■ \7g(p f ,R; e n ) = , (5) 

first derived by Eilenberger [^6| by eliminating the high-energy, short-distance structure of the full Green's function 
in Gorkov's equations. 27 The transport equation is supplemented by the normalization condition, 

g(p f ,R;e n ) 2 - -tt 2 1 , (6) 

which eliminates many unphysical solutions from the general set of solutions to the transport equation. |p6f 

The self energy, <t, has an expansion (Fig. 1) in terms of g (solid lines) and renormalized vertices describing the 
interactions between quasiparticles, phonons (wiggly lines), impurities and external fields. An essential feature of Fermi 
liquid theory is that this expansion is based on a set of small expansion parameters, small ~ ksTc/Ef , h/pf^o j ■■■ 1, 
which are the relevant low-energy (e.g. pairing energy) or long-wavelength (e.g. coherence length) scales compared 
to the characteristic high-energy (e.g. Fermi energy) or short-wavelength (e.g. Fermi wavelength) scales. p3 The 
leading order contributions to the self energy are represented in Fig.l. Diagram (la) is the zeroth-order in small and 
represents the band-structure potential of the quasiparticles. This term is included as Fermi-surface data for p/,Vf 
and Nf, which is taken from experiment or defined by a model for the band-structure. Diagram (lb) is first-order in 
small and represents Landau's Fermi-liquid interactions (diagonal in particle- hole space), and the electronic pairing 
interactions (off-diagonal in particle-hole space), or mean-field pairing self-energy (also the 'order parameter' or 'gap 
function'). Diagram (lb') represents the leading-order phonon contribution to the electronic self energy (diagonal) 
and pairing self-energy (off-diagonal); however, we confine our discussion to electronically driven superconductivity 
with a frequency-independent interaction. 

small 

a 



o 



small 1 



b b' 

9 9 9 



)— o— < 



dl d2 d3 

Fig. 1 Leading order contributions to the quasiclassical self-energy. 

In the spin-singlet channel, the order parameter satisfies the gap equation, 

A(p f ,R) = [ dp} V(p f ,p' f ) Tj2f(Pf,R\ e«) , (7) 

where f(pf,R;e n ) is the spin-singlet pairing amplitude and V(pf,pj) represents the electronic pairing interaction; 
this function may be expanded in basis functions for the irreducible representations of the point group, 

irrep d a 

v(p f ,Pf) = e v « £ y°* &) y ™ ■ ( g ) 
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where the parameter, V a , is the pairing interaction in the channel labeled by the ath irreducible representation, and 
the corresponding basis functions, {y a i(pf)\i — 1, ■ ■■,d a ], are orthonormal, (yai{pf)ypj{Pf)j — SapSij, where the 

Fermi surface average is defined by < A(pf) >f f = J dpf A(pf). 

Impurity scattering 

The summation of diagrams (Id) g ives the leading-order self energy from a random distribution of impurities in 
terms of the impurity t-matrix, |2q , |29| ] 

0j mj >(P/;en) = n imp i(pf,pf,e n ) , (9) 

i(Pf,Pf;tn) = u(ff,Pf) + N f J dp'} u{pfj"f) g(f'};e n ) t{p'},p};e n ) . (10) 

The first term is the matrix element of the impurity potential between quasiparticles at points pf and pf on the 
Fermi surface, rii mv is the impurity concentration, and the intermediate states are defined by the self-consistently 
determined quasiclassical propagator. 

For a spin-singlet superconductor with non-magnetic impurities, u(j>f,p}) = u(pf,pj)l, and the terms in &i mp 
that contribute in the transport equation lead to a renormalization of the Matsubara frequency and gap function; 
ie n = ie n — <Jimp{Pf] £n) and A(p/;e„) = A(j)f) + Ai mp (pf,e n ). Thus, the solution to the transport equation and 
normalization condition for the propagator becomes, 

g(p/;g») = ~ 7r ; = ^ ^=- (H) 



\Jt n {pfi£n) 2 + |A(p>;e„)| 5 



In the second-order Born approximation for the impurity t-matrix (this is not essential, but simplifies the following 
discussion) . The impurity renormalization of the off-diagonal self-energy is given by, 

A(p);e„) = A(p» + / dp[w(p f ,p}) . , (12) 



where w(pf,pj) = 2irni mp Nf\u(pf ,pf)\ 2 is the scattering rate in the Born approximation. Note that the integral 
equation for A(pf;e n ) has the mean-field order parameter, A(pV), as the driving term. The scattering rate w(pf,pj) 
has the full symmetry of the normal metal; thus, it too can be expanded in basis functions for the irreducible 
representations of the point group, 



wip f ,p'f) = E E y«i&f)y*<?}) > ( l3 ) 



a: 



where l/2r Q is the scattering rate for channel a. The integral equation for the renormalized order parameter separates 
into algebraic equations for each representation, 

a.= / dp f y* ai (? f) A(? f)+ ± [ d p f ^MMmM . d4) 



The driving term is non-zero only for the irreducible representation corresponding to A(p/). Thus, the resulting 
solution for the impurity renormalized order parameter necessarily has the same orbital symmetry as the mean-field 
order parameter; and the magnitude of the impurity renormalization is determined by the scattering probability for 
scattering in the same channel as that of A(pf). |3(j] The argument also holds for the full t-matrix. 

For isotropic ('s-wave') impurity scattering the renormalized Matsubara frequency and order parameter become 



e« = e« + — ( , = ) , (15) 

' !<? n + \AW;e n )\2/ 

V'f 
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A(p / ;e„) = A (p -» + i-( ) . (16) 

2r \^ + |A(p/;e„)p/^ 

Thus, for an s-wave order parameter these equations give identical renormalization factors for both the Matsubara 
frequency and the order parameter, i.e. 

e n A(e ra ) 1 1 , > n7 , 

in which case the impurity renormalization drops out of the equilibrium propagator and gap equation. |?l|]32] | However, 
s-wave superconductors are exceptional; for any unconventional order parameter impurity scattering is pairbreaking. 

E 

Consider an unconventional superconductor with impurities in which the scattering is dominated by the identity 
representation. If there is an element of the point group, 1Z, which changes the sign of A(pj), i.e. A(p/) — > — A(p/), 



then from eq.([16|) the impurity renormalization of the order parameter vanishes identically: A(jf/; e„) = A(p/). The 
cancellation between the impurity renormalization factors for the Matsubara frequency and order parameter no longer 
occurs, with the consequence that impurity scattering suppresses both T c and the magnitude of the order parameter. 

For isotropic impurity scattering (not restricted to the Born approximation), the renormalization factor for the 
Matsubara frequency, e„/e„ = Z(e n ), is independent of position on the Fermi surface and given by 

7U \ - 1 +F Z(e n )V(e n ) 

ctn z (d ) + (Z(e n )e n V(e n )y 

with 

(19) 

where r„ = n imp /irNf and S Q — tan _1 (7riVyito) is the s-wave scattering phase shift in the normal state. In the Born 
limit, 5q — » TrNfUQ, while in the strong scattering limit (Nfiio — > oo) we obtain the unitarity limit, 5q — > 7r/2. 

Given the gap function, A(p/), the impurity renormalization is easily calculated. The magnitude and temperature 
dependence of the order parameter are calculated self-consistcntly from the mean-field gap equation, 

A(p f ) = f dp' f V(p f ,p}) nT ' £ ' .ffiL^n,. ■ (2°) 

The linearized gap equation determines T c in terms of the pairing interaction, frequency cutoff lo c and impurity 
scattering rate. At T c only the dominant pairing channel a is relevant and the linearized gap equation becomes, 



|e n |<u 



1 



v a ^ |e„|+r' 

where r = r„ sin 2 <5o is the pair-breaking parameter, 

( T u sin 2 S — %f = -Kn imp N f ul , (Born limit) 



(21) 



~~ \ r. u = , (unitarity limit) . 1 ' 

For r = 0, this equation determines the clean-limit value of the transition temperature, T co . Eliminating the pairing 
interaction and cutoff gives the well-known Abrikosov-Gorkov formula, []33f except that the pairbreaking parameter 
is determined by non-magnetic scattering, |6| 

*a+^)-*m=i»o). (23) 

where tJj(z) is the digamma function. 



2 27rT r J V 2 J \T, 
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Finally, the linearized gap equation is used to eliminate the pairing interaction and cutoff in favor of T c in the full gap 
equation (|20|). For pairing in a ID representation, V(pf,Pf) — V e{pf)e{pf)* , the order parameter is A(p/) = Ae(p/), 
where A is obtained from eq. (p0|) . Multiplying eq.(20) by e(p/)*, integrating over the Fermi surface, and adding and 
subtracting the RHS of the linearized gap equation (with T c — > T) to eliminate V gives, 



ln(T/T )+^(i 



r 1 r 

2itT j n 2 2]rT/ 



|e(p/)| 2 



pi 



(24) 



which is solved self-consistently with eq.(|l§|) to give Z(e„) and A as a function of T/T c and T. 



Gauge-invariant coupling to the condensate flow field 



The self-energy term representing the diamagnetic coupling of quasiparticles to a static magnetic field (Fig. lc) is 
determined by local gauge invariance. Under a gauge transformation of the Fermion fields, ip(r) — > ijj{r)e~ l ^ r ^ 2 , the 
quasiclassical propagator transforms as 



g^g' = #(A)t g <7(A) 



(25) 



where U(A) = exp(+|A(i?)f3), as does the self-energy and order parameter. |2j| Applying this transformation to 
the transport equation (||) gives ?/(A)^ Q[g, a] U{A) — Q[g' ,<j' + Thus, the form of the transport equation is 
invariant, but the local gauge field generates an additional self-energy, dy = — iU(Ay Vf(f>f) ■ V U(A). This property 
of the transport equation is used to eliminate the phase degree of freedom of the order parameter in favor of a spatially 
varying flow field. We parametrize the spatial variations in terms of a physical gauge, the phase x(^)i an d the local 
amplitude, A (p f ;R) = \A(R)\ e(p» 



A(f f ,R) = U[x(R)} A Q (p f ,R) U^x(R)} ■ 
Thus, the transport equation becomes 

ie n T 3 - A - a v - a , g' + ivf(pf) ■ Vg = , 



(26) 
(27) 



where a v = kvf(Pf) ' ^X^3- The diamagnetic coupling to a magnetic field, b = V x A, is then determined by gauge 
invariance, and can be represented in terms of the gauge-invariant condensate flow field, 



;(Vx 



2e 



A). 



and the self-energy, 



VfiPf) ■V s (R)t 3 



(28) 



(29) 



In the Meissner geometry (H parallel to the interface) the driving term associated with the applied surface field is 
of order, 



irT c 



H 
K 



(30) 



where H c ~ 0o/(£^) is th e thermodynamic critical field. This term should be compared to the gradient term arising 
from spatial variations of the screening current, 



v f ■ Vg 




v f ■ V(ct„/A) 


ttT c 




7rT c 



LE. 



(31) 



In the strong type II limit the velocity field is effectively uniform on the scale of the coherence length so we are 
generally justified in dropping the gradient term in eq. (27). 
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Linear response 



The linear response limit is simply obtained from a perturbation expansion of the propagator, transport equation 
and normalization condition, and is expected to be valid for low magnetic fields, \<j v /t:T c \ ~ H / H c -C 1. Assume an 
expansion of the form, g = go + g\ + where go is the zero- field solution to the transport equation given by eq.(^Tj), 
and gi is the first-order correction to the propagator, formally of order \gi\ ~ 0\(a v / A)go\. The linearized transport 
equation and normalization condition, 

ie„f 3 -A,<?i -[a v ,g o ] = , {go,9i} = 0, (32) 

are inverted with the aid of eqs. (O) and (0) to give 

. |A| 2 - zg»f 3 A 
gi = na v = — . (33) 

The resulting supercurrent calculated from eq. (@) can be written in terms of the superfluid density tensor, 

(p s )^ = 2N f f dp f $(p» v){p f ) v](p f ) , (34) 



where 

which reduces to the angle-dependent Yosida function in the clean limit. 

Fermi-liquid effects 

Fermi-liquid effects arise from the leading order electronic self-energy (diagram lb). For the diamagnetic response 
the most important Fermi-liquid effect is the screening correction to the diamagnetic current; the relevant self energy 
is, a b (pf,R) = abipf, R) f 3 , with 

a b (p f ,R) = f dp} A™ r (p f ,p})Tj2 9(Pf, R; e n ) , (36) 

where A cur (pf,pf) is the dimensionless quasiparticlc interaction. 

Fermi liquid effects can contribute substantial temperature-dependent corrections to the penetration depth as the 
gap opens and the number of thermal quasiparticles drops rapidly. We include formulas for the Fermi-liquid correction 
to the supercurrent for a model uniaxial Fermi surface. The position on the Fermi surface can be parametrized by 
the direction of the Fermi wavevector, p, and the Fermi velocity is given by, 

Vf = v^J (p x x + p y y) + vfp z z . (37) 

Similarly, the quasiparticle interaction is parametrized by two Landau parameters corresponding to current flow in 
the basal plane and along the z-axis, 

A cur (Pf,p}) = AH (p x p' x + p y p' y ) + A 1 - {p z p' z ) . (38) 



The linear response result for the supercurrent is easily obtained from eq.(|33j) with the replacement a v —> a sc — 
Vf(pf) ■ v s + at,(j)f). The resulting current is given by 



MR) = -2eN f j dp f v f (p f ) $(p» a sc (Pf) , (39) 
with the screening field satisfying the self-consistency equation, 
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°sc{Pf) = Vf{p f ) -v s + J dp} A cur {p f ,p}) a sc (pf) . (40) 

Equations ( p4p^ , |3^ , ^0| , p^| , p4[ ) are the basic equations used to calculate the temperature-dependent penetration depth 
for unconventional superconductors. 

For a gap |A(p/)| with four- fold rotational symmetry about the z-axis the resulting supercurrent is given by a 

diagonal superfluid density tensor with in-plane (pi') and z-axis (pjr) superfluid densities given by 

^^ ■■iiV ' (41) 

with 

h = lf§ (pi+Pi) m), ^=3 J §Ptm)- (42) 

Since <p\\,± ~ Aq ~ (1 — T/T c ) near T c , the Fermi liquid renormalizations of the penetration depth drop out near T c , 
but may give substantial corrections to pi' ,_L (T) at low-temperatures. 34 3^j 



III. ZERO-FIELD PENETRATION DEPTH OF A D X 2_ 



SUPERCONDUCTOR 



Consider the model of the CuO superconductors based on a d x 2_ y 2 order parameter. The general form of the order 
parameter is, 



A(Pf) = A{pl-p 2 y )*I(p f ), 



(43) 



where I(pf) is invariant under the full point group, and p X)V define the direction of the Fermi wavevector in the basal 
plane of the Fermi surface. Note that the nodes are required by the broken reflection symmetries. Two parameters 
determine the excitation spectrum in the clean limit; (i) the maximum value of |A(p/)| (= Ao), and (ii) the angular 
slope of the gap near the node, 



1 d|A(tf)| 



A dtf 



(44) 



where i? is the angle measured relative to one of the nodes (see Fig. 2). 
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Fig. 2 Gap function for a d x 2_ y 2 superconductor. 







A simple two- parameter model for |A(p/)|, which is useful for numerical calculations, is 

The maximum gap, for a fixed fj,, is obtained from a self-consistent solution to the gap equation. Figure 3 shows 
solutions of eq.(|24|) for Ao(T) as a function of temperature and impurity scattering. We obtain a gap ratio of 
Ao/T c = 1.9 at T = for /j, = 2.7 and r = 0. Note that the leading temperature-dependent correction to the gap 
parameter for T <C T c is JAq(T) ~ T 3 , in the clean limit (appendix A). 
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Fig. 3 Maximum gap for a d x 2_ y 2 superconductor as a function of temperature. The inset shows the T deviation 
of Ao(T) at low-temperature. 



The maximum gap is relatively insensitive to the angular slope of the gap near the nodes, except for small /i, 
in which case the nodal region occupies a significant fraction of phase space. This behavior can be qualitatively 
understood by noting that the condensation energy at zero temperature is effectively determined by the strength of 
the pairing interaction, and therefore T c . For an anisotropic gap A is enhanced to compensate for the regions of 
small gap. If we assume that the Fermi-surface average of |A(z9)| 2 is constant (fixed by T c ), then the maximum gap 
is given by Aq(/i) ~ 1.8T C /(1 — 8/37r/i), which is qualitatively the behavior obtained from the numerical solution to 
the gap equation shown in Fig. 4. Note that the ususal one-parameter d x 2_ y 2 gap, A = A (|5 2 — p^), corresponds to 
[1 = 2. 




0.00 2.00 4.00 6.00 8.00 10.00 12.00 



p. 

Fig. 4 Maximum gap for a d x 2_ y 2 superconductor as a function of the angular slope parameter, fi. 
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The suppression of T c by impurity scattering is determined by the pair-breaking parameter T = T u sin 2 So- Equation 
( ^3| ) for T c implies that unconventional pairing can be sensitive to impurity scattering, e.g. the superconducting 
transition is completely suppressed for (r/7rT co ) crit = — 0.28, which corresponds to an impurity mean-free 

path of l cr it — Vf/2T ~ 3.6^0j where £o = Vf/2irT c o is the coherence length in the clean limit. The relatively small 
coherence length in the CuO superconductors is then an advantage for an unconventional order parameter. For weak 
pair-breaking, T/2t:T c <C 1, the suppression of T c is given by, AT c /T c q = —ttT/8T c q. Thus, for CuO superconductors 
with T c o = 100-ftT and a suppression of less than Q.5K we have T < 1.3K. For an in-plane coherence length of £o = 14A 
this corresponds to an impurity mean-free path, / > 3, 450A.Q 

The magnitude of the gap parameter is also suppressed by impurity scattering. Figure 5 shows the suppression of 
Ao(0) as a function of T. Note that Aq(0), in contrast to T c , is more strongly suppressed in the unitarity limit than 
in the Born limit for the same suppression of T c . 




Fig. 5 Maximum gap at T = 0.01T C for a d x 2_ y 2 superconductor as a function of the impurity pair-breaking 
parameter. The inset shows the impurity pair-breaking effect on T c . 



In the clean limit the angle-resolved density of states, obtained from eq. (|Tl]), is given by the familiar BCS form, 

Af(p f ,E) = lE [ G(E* - |A(p»| 2 )- (46) 

V E - \ A (pf)r 

For low energies, \E\ <C Ao, the total density of states is dominated by the low-lying states near the nodes, and is 
linear in \E\. For the model gap function in eq.(|45|) the total density of states is 

ST{E) = J dp f N(p f ,E) ~ ffl , \E\<A . (47) 

These low-lying excitations are responsible for the linear temperature dependence of the penetration depth for a 
superconductor with a line of nodes (or point nodes in 2D). Not surprisingly, the density of low-lying states near 
the nodes is determined by the angular slope of the gap; therefore, \x also determines the coefficient of the linear 
temperature dependence of the penetration depth (in the clean limit). 

Although the small pair size in the cuprates leads to relatively weak suppression of T c from impurity scattering, 
the density of states at low-energy, \E\ <C Ao, and therefore the leading temperature dependence of the penetration 
depth, are more sensitive to impurity scattering, particularly in the strong scattering limit. Figure 6 shows the density 
of states as a function of the scattering phase shift for fixed T u = O.IAq, corresponding to a rather high impurity 



2 For non-magnetic, s-wave impurities in 2D, the impurity resistivity is given by p~^ p = e 2 Nfv'j/AF, where F is the same pair- 
breaking parameter that enters the Abrikosov-Gorkov formula. Thus, we can express T c /T c q as a function of pi mp , independent 
of the scattering phase shift, do- However, inelastic scattering, which is important at T ~ T c , destroys this simple result. [§6l|37| 



11 



concentration. The modification of the density of states at low-energy is negligible in the Born limit (So — 7r/20) and 
remains essentially linear even for intermediate phase shifts (So = tt/4). However, as the strength of the scattering 
increases a finite density of states at E = develops, becoming of order Af(0) ~ 0.4 in the unitarity limit. In addition, 
Af(E) deviates from linearity below E ~ 0.4Ao. This is the cross-over energy scale, below which impurity scattering 
strongly modifies the low-energy spectrum. The cross-over energy, e* can be calculated from the lowest energy scale 
for the renormalized Matsubara frequency e n as T — + 0. From eq. ( |l8|) the cross-over scale at T = is given by 

<tT u (l-a)+a(e*V(e*))*' { ' 

where a = sin 2 So, and T>(e*) can be calculated to leading logarithmic accuracy (ln(Ao/e*) ^> 1) for the d x i_ y i state, 

P( £ *)^^-ln(2A /£*). (49) 



In the unitarity limit this scale can be a sizable fraction of Ao even in the dilute limit, e* ~ ^/T\A~o, but in the 
Born limit the cross-over scale is exponentially small, e* ~ 2Ao exp(— /i7rrAo/2). For more detailed discussions of 
the density of states of d-wave superconductors see Refs. ( p8|,^9|). 




E/A 



Fig. 6 Density of states vs. scattering phase shift. The impurity concentration is fixed with r„/Ao = 0.1. Note 
the finite density of states and the suppression of the maximum gap in the unitarity limit. 



Penetration depth for the d x 2_ y 2 gap 

The temperature dependence of the penetration depth for a d x i_ y i gap function is obtained from eq.(|34|) (neglecting 
Fermi liquid corrections) . We assume a cylindrical Fermi surface and a gap function parametrized by /i and Ao as in 
eq. (p5|) . The in-plane penetration depth becomes 

1 _^e 2 N f (v\) 2 rdd |A(0)| 



An 



where the integration is over the Fermi circle in the basal plane, and Z(e n ) is the impurity renormalization factor for 
s-wave scattering centers. In the clean limit the leading correction to the penetration depth at low-temperatures is 
linear in T, typical of a gap with a line of nodes, J3^ , fifi| 442|| 



5A|i(T) / 2 In 2 



Aii (0) \d\A(<&)\/d#\ nodt 



T+... ■ T«T C , (51) 



where d\/S.(d)\/ d§\ no d e is the angular slope of |A(i?)| at the node. For the model gap in eq. (|45|), d|A(i?)|/<ii?| no( j 
fj,A . 
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The report by Hardy, et al. Q| of a linear temperature dependence to <5A||(T) for single crystals of YBa2Cu30g.95 
provided the first substantial evidence for a superconducting state with a line of nodes in the excitation gap. Figure 
7 shows the low-temperature data for <5A|| (T) reported in Ref.( ff3|). The solid line is a calculation of the penetration 
depth for a d x 2_ y 2 gap with the angular slope adjusted to fit the data for T < 20 K. For the absolute penetration 
depth we assume A|| (0) = 1, 400 A. (ij] The slope of the penetration depth at low temperature is 4.3 A/K, and the fit 
to eq. (pl|) gives [i = 2.7, A o (0)/T c = 1.9, and d|A($)|/d$| no( 2 e = /iA o (0) ~ 5.1 T c . For comparison, the one-parameter 
d x 2_ y 2 model, A = Aq cos(2</>) with Aq calculated self-consistently, gives d\/S.( , d)\/ d-d\ no d e = AiAo(O) = 4.3 T c . 




Fig. 7 Penetration depth as a function of temperature. The (*) are experimental data obtained from Ref. ( 
0) the solid line is the best fit with /i ~ 2.7, Aq/T c ~ 1.9. The inset shows the effect of a tiny gap at the nodal 
positions on the temperature dependence of the penetration depth. The axes are the same as those of the main 
graph. The middle (top) curve in the inset corresponds to a tiny gap of A m i„/Ao = 1% (2%). 



Also shown in Fig. 7 is the effect of a tiny gap at the nodal positions. Tiny gaps can arise in strongly 
anisotropic conventional superconductors, |45|] or if the unconventional order parameter contains a small compo- 
nent of another representation. For example, an order parameter of the form d x 2_ y 2 + i\fed xy has a gap function, 
|A(p/)| = Ao[cos 2 (2</>) + e sin 2 (2(f>)] 1 ^ 2 , which is strongly anisotropic for small e with a tiny gap of order eAo near the 
nodes of the d x 2_ y 2 component. [Q The data of Ref. ( Q|) implies that e < 2%. 

Precision measurements of the penetration depth in thin films show a T 2 behavior at low temperatures rather than 
the linear temperature dependence characteristic of a d x 2_ y 2 order parameter. [^7[- p2[p^ ] This difference may be due 
to scattering by a higher concentration of defects present in the films. Hirschfeld and Goldenfeld |4|] argue that the T 2 
dependence in films and the T dependence of <5Aii(T) reported for single crystals can be understood within the same 
d x 2_ y 2 model for the pairing state provided the films are relatively dirty compared to the single crystals. However, 
impurity scattering or defect scattering is pair-breaking in unconventional superconductors, so in order to explain the 
weak or negligible suppression of T c in films (compared to T c in the single crystals of Ref.( ji]), the authors of Ref. 
( |l|]) argue that the scattering responsible for the T 2 dependence of <5A||(T) results from a dilute concentration of 
strong scattering centers with phase shifts near the unitarity limit. A small concentration of unitarity scatterers leads 
to strong modification of the density of states in the small phase space region near the nodes. Since <5A|| is determined 
by these low-energy excitations, the temperature dependence of 8\\\(T) for T <C T c may be strongly modified even 
when the suppression of T c by a dilute concentration of scatterers is negligible. 

In the Born limit (weak scattering) the density of states, even near the nodes is nearly unchanged. Thus, a much 
higher concentration of defects is needed to generate <5A|| (T) ~ T 2 for T < 0.2T C , which is accompanied by a sizeable 
suppression of T c . The sensitivity of the low-energy excitation spectrum to the scattering strength is reflected in the 
temperature dependence of the penetration depth shown in Fig. 8. The concentration of s-wave scatterers is fixed 
and the curves show the evolution from a linear T dependence in the Born limit (So = 7r/20), and intermediate phase 
shifts, to the T 2 dependence in the unitarity limit (5q = tt/2). Note that the cross-over from SXu ~ T to #Am ~ T 2 is 
abrupt, occuring very near the unitarity limit for dilute concentrations. A sharp cross-over in the excitation spectrum 
as a function of phase shift has also been noted by Preosti, et al. pq |. Thus, it is worth emphasizing that for dilute 
point impurities it is not merely 'strong scattering' that is required to obtain <JAm ~ T 2 with minimal reduction in T c , 
but scattering with So — ► 7r/2. 
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Fig. 8 Penetration depth as a function of temperature for phase shifts ranging from the Born limit to the 
unitarity limit. The concentration is fixed with F u — 0.01 T c . 



In the unitarity limit even dilute concentrations of impurities can strongly modify the low-energy spectrum. In 
the unitarity limit the cross-over energy scale is e* ~ y/T u Ao. Thus, SXn(T) deviates from the linear T dependence 
for T < e* . The sensitivity of SXu(T) at T <C T c for unitarity scattering (see Fig. 9) places a strong constraint 
on the concentration of scatterers that can be present in clean single crystals that shows SXn(T) ~ T down to low 
temperatures. For the data of Ref. ( 0) we find r u /A < 0.0001 given that SX\\(T) ~ T down to T ~ 1 K. 




Fig. 9 Penetration depth as a function of temperature and impurity concentration (~ r«) in the unitarity 
scattering limit. 



IV. NONLINEAR CURRENT RESPONSE 



In the Meissner geometry the screening current is proportional to the applied surface field, j s ~ cH/ A. As H is 
increased nonlinear field corrections to the constitutive equation for the supercurrent may become significant. In 
conventional type II superconductors nonlinear corrections to the current-velocity relation arise from the thermal 
population of quasiparticles, and vortex nucleation generally occurs before these nonlinear effects become important. 

In unconventional superconductors with nodes in the excitation gap the nonlinear field correction to the supercurrent 
is substantially larger than in conventional superconductors with other similar material properties. The origin of the 
anomalous nonlinear Meissner effect is the contribution to the screening current associated with the quasiparticle 
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states near the nodal lines. As a result the nonlinear Meissner effect may be used to detect the nodal structure of the 
gap of an unconventional superconductor. |l5f| 

In the limit \vf ■ V|A||/7rT c ~ £/A <C 1 the current can be expressed as a local function of the condensate velocity, 
v s (R). In the presence of a condensate flow the local solution to the transport equation, g(pf,R; e„), is given by eq. 
(pi]), but evaluated with ie n — * ie n — a v , where a v (pf,R) — Vf ■ v s (R). The current response obtained from eq. (^) is 

I = ~2eN f f dp f v f (Pf) *T£ °vtff)-iin (52) 

J ^(e n + ia v (p f )y + \A(p f ;e n )\2 

One point to note in calculating the density of states, or current, at finite flow with impurity scattering is that the 
impurity-renormalization of A(p/) need not vanish, even for s-wave impurity scattering and (A(p/)) = 0. The reason 
is clear from eq. (|l^) for A(p/ ; e„); the kernel no longer vanishes by symmetry with the replacement, e n — > e n + ivf-v s 
for a general flow field. However, for special directions of v s , e.g. v s parallel to a node, the impurity correction to 
A(p/) vanishes . A = A(pf). 

Equation ( p2|) for the current can be transformed by contour integration and analytic continuation to the real axis 
to give 

r r+oo 

I = -2eN f / dp f v f {p f ) / dE f(E) W+(p f ,E) - M. (p f ,E)} , (53) 

where M±{pf , E) is the density of states for quasiparticles that are co- moving (+«/ • v s > 0) and counter-moving 
{—Vf ■ v s < 0) relative to the condensate flow. The integral is taken over the half space a v = Vf ■ v s > with the 
counter-moving excitations included by inversion symmetry: Ci<j v = —o v and |A(Cjp/)| = |A(p/)|. This result is 
general enough to cover nonlinear field corrections to the current for superconductors with an unconventional order 
parameter and pair-breaking effects from impurity scattering. 




A|-V f -V s -|A| -|A| + V f -V s |A|— V f -V s |A| A|+V f -V s 



Fig. 10 Density of states for co-moving (+Vf ■ v s ) and counter-moving (— Vf ■ v s ) excitations at a point pf on the 
Fermi-surface where \vf ■ v s \ < |A(p/)|. 



The difference in the nonlinear current-velocity relation for conventional and unconventional order parameters 
appears in the contributions to the current from the co- moving and counter- moving excitation spectrum at T = 0. 
The spectrum is shown Fig. 10 in the clean limit for a specific direction pf in which Vf ■ v s < \A(pf)\. At zero 
temperature only the co- moving and counter- moving quasiparticle states with E < contribute to the current. 



Nonlinear current: conventional gap 

For a conventional superconductor with an isotropic gap at T = the current is easily calculated from the difference 
in the number of co-moving versus counter-moving quasiparticles that make up the condensate, 

j s = -2eN f / dpfVf(p f ) [2vf ■ v s ] = ~epv s , v s < v f /A , (54) 

Jcr„>0 

with p = Nfv 2 j for a cylindrical Fermi surface. The main point is that the current is linear in v s for velocities up 
to the bulk critical velocity, v c — Vf/A^. At v s — v c the edge of the spectrum for the upper branch {E > for 
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v s < v c ) of counter-moving excitations drops below E = 0, and the edge of the spectrum for the lower branch (E < 
for v s < v c ) of co-moving excitations shifts above E = 0. As a result the current carried by the condensate drops 
rapidly above the critical velocity. The current is nonanalytic at v c because a branch of counter-moving (co-moving) 
excitations that are unoccupied (occupied) for v s < v c become occupied (unoccupied) v s > v c . For example, for a ID 

Fermi surface the current becomes, j s = —2eNjvj |u s — Q(v s — v c ) y/v^ — v%\. 

At non-zero temperatures thermal occupation of the upper branches and de-population of the lower branches reduces 



the condensate supercurrent. In the clean limit eq.(53) can be transformed to 



2eN f / dp f v f (v f -v s )+j qp , (55) 



Jqp 



= -4eNr 



f r <% f dp f v f [f (Ve+W + v f ■ v s ) 



(56) 



which separates the condensate contribution to the current, the fully occupied negative energy branches shown in 
Fig. 10, from j qp , the current carried by the excitations associated with population of the upper branches and de- 
population of the lower branches in Fig. 10. The current carried by the excitations is a backflow current. For low 
velocities the net current is linear in v s , j s — ~ ep s (T) v s , where p s {T) is the superfluid density in the two-fluid model, 



p s {T) = p-Nfv) 



sech 2 ^ 2 + jAj 
2T 



(57) 



For larger flow velocities the linear relation breaks down. The leading correction to the current-velocity relation for 

v s < v c is, 



js = -ep s (T) v s 



1-aCO - 

V c 



(58) 



The nonlinear correction is third-order in v s and determined by the coefficient a(T) > and the bulk critical 
velocity, v c — A(T)/vf. There are two sources to the nonlinear current. At finite v s there is a difference in the 
thermal occupations of co-moving and counter-moving excitations. Since the counter-moving branch of has a larger 
occupation than the co-moving branch (f(E — VfV s ) compared to f(E + VfV s )), the thermal excitations further reduce 
the current compared with the linear response value. In addition the condensate velocity is also pairbreaking, i.e. v s 
reduces the magnitude of the mean-field gap parameter, further reducing the current density at finite flow. Figure 11 
shows the temperature dependence of the nonlinear correction to eq.([38|). Note that a(T) ~ exp(— A/T) for T — > 0. 
There are no excitations that contribute to the backflow current at T = and the constitutive equation is strictly 
linear for velocities below the bulk critical velocity, v c = A/uy. 
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Fig. 11 Temperature dependence of p s {T) and the nonlinear coefficient a(T) for an s-wave gap. 
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The relevance of the nonlinear current-velocity constitutive equation to the penetration of magnetic fields into a 
superconductor is qualitatively clear. The reduction of the current by the flow reduces the effective superfluid density, 
and, therefore, increases the penetration of the field into the superconductor. Since the current is proportional to the 
field in linear order, the leading correction to the effective penetration length is quadratic in the surface field. Solution 
of the nonlinear London equation in the Meissner geometry gives the following result for the field-dependence of the 
penetration depth (see appendix B), 

-J— = J- !i- 3 - 

X(T,H) X(T) \ 4 

where 1/A(T) 2 = 87re 2 p s (T)/3c 2 is the zero-field London penetration depth, and Hq(T) — e \(T)/c v c (T) is of order 
the thermodynamic critical field. Thus, in the London limit nonlinear Meissner effect in a conventional superconductor 
is exponentially small at low temperatures and is quadratic in H / Hq. 



H 



Ho(T) 



(59) 



Nonanalytic supercurrents at T=0: d x 



gap 



The expansion of the current in v s breaks down for an unconventional superconductor with nodes in the gap. This 
is clear from eq. (52); a Taylor expansion in a v fails for T <C T c when there are directions pf where |A(p/)| and |e„| 
are always small compared to |cr„(p/)|. In the clean limit for a d x 2_ y 2 gap the breakdown of the Taylor expansion 

leads to a nonanalytic current-velocity relation at T = of the form j s = —epv s {1 — |w s |/uo}, where vo ~ Ao/u/.f^] 
The physical origin of this nonanalytic current is easily understood by considering the d x i_ y i gap with v s directed 
along a node, as shown in Fig. 12 (left panel). 





Fig. 12 Phase space contributing to the quasiparticle backflow jets at T — for v s \\ node, and v s \\ antinode. 



For any v s ^ there is a region of the Fermi surface with |A(p/)| + Vf ■ v s < 0, in which the upper branch of the 
counter-moving excitations (see Fig. 10) have negative energy and become populated. The non-analytic dependence 
on v s reflects the occupation of this counter-moving branch of excitations at T = 0. Figure 12 illustrates the phase 
space contributing to the backflow current. For \vfV s \ <C Ao the wedge of occupied states is — i? c < ■& < i? c , with 
•d c = v*fV s /nAo, where v*j is the Fermi velocity at the node and /xAo is the angular slope of |A(#)| at the node. The 
current is calculated by transforming cq. (153) to 



-2eN f 



dp f vf 



(v f -v s ) + 2 dU(E(0+v r v s ) 



(60) 



3 Nonanalytic currents have been investigated in superfluid 3 He-A, initially by Volovik and Mineev. This work is closely 
related to a number of theoretical investigations of the hydrodynamical equations of superfluid 3 He in the limit T — > (see 
Ref. ( J54|]). For an analysis of the non-analytic current in 3 He-A see Ref. ( |5q| ); these authors also calculate the non-analytic 
current- velocity relation for an axially symmetric, polar state with A ~ p z . The polar model is examined in more detail by 
Choi and Muzikar. Q 
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where E(£) = \/£ 2 + |A(p/)| 2 . The first term is the condensate current, —epv s . The backflow current at T — is 
easily calculated from the phase space of occupied counter- moving excitations. With the velocity directed along the 
nodal line p x = p y , i.e. v s = v s x' as shown in Fig. 12 (left panel), the occupied states give 



j qp = -2eN f [dp f Vf <d(-v f ■ v s - \A(p f )\) yj (v f ■ v s )2 - \ A(p f )\< 



(61) 



d() 



-•IrVjV) J ~— yJ(v* f V 3 )* ~ ( M A ^) 2 (-£') . 



To leading order in (v%v s /Aq) we obtain a total current of 



-epv s 



1 - 



(v s 1 1 node) 



(62) 



where vq = /iAq /v% is of order the bulk critcial velocity scale. Equation (|62|) clearly holds for v s directed along any 
of the four nodes. Note that the current is parallel to the velocity, and that the counter-moving excitations reduce 
the supercurrent, as expected. Also the nonlinear correction is quadratic rather than cubic, as is obtained for the 
conventional gap, and with the characteristic scale determined by vq — fxAo/v%. 

Unlike the linear response current, the nonlinear quasiparticle current is anisotropic in the basal plane. A velocity 
field directed along the maximum direction of the gap (antinode), v s — v s x, produces two counter-moving jets, albeit 
with reduced magnitude because the projection of v s along the nodal lines is reduced by l/y/2 (Fig. 12, right panel). 
The critical angle defining the occupied states in this case is given by -^(v s /vq). The resulting current is easily 
calculated to be 



Jqp 



-ep 



v s /V2 



v 



MO + 



v 



giving a total current of 



Js 



-ep v s 



1 \= \ (v s 1 1 antinode) . 

V2 vo 



(63) 



(64) 



which is again parallel to the velocity and has a quadratic nonlinear correction. However, the magnitude of the 
nonlinear term is reduced by l/\2. This anisotropy is due to the relative positions of the nodal lines and is insensitive 
to the details of the anisotropy of the the Fermi surface or Fermi velocity because the quasiparticle states that 
contribute to the current, for either orientation of the velocity, are located in a narrow wedge, i? < i? c ~ (v s /v ) <C 1, 
near the nodal lines. Thus, the occupied quasiparticle states near any of the nodes have essentially the same Fermi 
velocity and density of states; only the relative occupation of the states is modified by changing the direction of the 
velocity. 

The dependence of the supercurrent on the positions of the nodal lines in momentum space suggests that the 
anisotropy can be used to distinguish different unconventional gaps with nodes located in different directions in 
momentum space. For example, the d xy state (B2 S representation), A ~ p x p y , would also exhibit a four- fold 
anisotropy, but the nodal lines are rotated by 7r/4 relative to those of the d x 2_ y 2 state. The order parameter A ~ 
Px Py{p x — Py), corresponding to the A 2g representation, would exhibit a more complicated 2 x 4— fold anisotropy. 

Anisotropy in the in-plane current implies a similar anisotropy in the field dependence of the in-plane penetration 
length. Consider the geometry in which the superconductor occupies the half-space z > 0, with z\\c. For a surface 
field H directed along a nodal line, Maxwell's equation combined with eq. (|6^) for the current and the gauge condition 
V • v s = reduces to 



47re 2 



dz 2 



Js[Vs 



=-i (1 - 



Vo 



(65) 



We define the effective penetration length in terms of the static surface impedance, 1/A|| = —(l/H)(db/dz)\ z =o- The 
solutions for both half-space and thin film geometries are discussed in appendix B. We obtain, 



11/ H 

A N (#) ~ h \ ~ 



if 1 1 node . 



(66) 
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and similarly for fields directed along an antinode, 



where Am is the zero- field penetration depth at T — and Hq = 3ci>o/2eA|| is the characteristic field scale. Using 
v = pA /v* f , l/X 2 = 4ne 2 N f v 2 f /c 2 , £ N = v f /2nT c and H 2 /8n = ±N f A 2 gives Hq ~ %p{v f /v* f ) 3 H c ~ 0(H C ), where 
H c is the thermodynamic critial field, vj is the rms average of the Fermi velocity and is the Fermi velocity at the 
node. We estimate H c (0) ~ 8.5 kG and Hq ~ 3.4 T if we neglect the anisotropy of the Fermi velocity. Fermi surface 
anistropy can change the characteristic field significantly. Using a next-nearest neighbor tight-binding model for the 
Fermi surface fit to the LDA result for YBa2Cu307_5, |37| we estimate the anisotropy to be v*^/vf ~ 1.1, and gives 
H ~ 2.5 T; however, the anisotropy of Vf is sensitive to the hole concentration near half filling. Assuming H = 2.5 T 
and H cl (0) ~ 250 G appropriate for twinned single crystals of YBCO, ||| then over the field range, < H < H cl , the 
change in A||(ff) is of order SX\\ ~ A N H cl /H ~ (1, 500 A)(250 G)/(2.5 x 10 4 G) ~ 15 A. The magnitude of 5\\\{H) 
also depends on the measurement technique. If 5A\\(H) is measured from the inductive response of a low-frequency 

a.c. field, A u , in the presence of a parallel d.c. field, v s , then from eq. ( |6^ ) Sj^ = — ^-p{l — 2 \v s \/vq} A^; and the 
change in the penetration depth with the static field is a factor of 2 larger than the d.c. result in eqs. ( |66| ) and (|6^). 

Observation of the anisotropy and linear field dependence of <5A||(-ff) could provide strong support for a d x 2_ y 2 
order parameter. Below we consider thermal and impurity effects which might mask or wash-out the characteristic 
anistropy and linear field dependence characteristic of pure material at T = 0. 

Thermal excitations: cross-over to analytic behavior 

At finite temperatures thermally excited quasiparticles occupy the upper branch of the counter-moving band shown 
in Fig. 10, and for T <C Aq these thermal quasiparticles are predominantly in the nodal regions. In the limit 
VfV s -C ttT the thermal excitations dominate and the nonlinear corrections can be obtained from a Taylor expansion 
in (vfVs/nT). In the opposite limit nT -C VfV s -C 7rT c , the non-thermal jets that give rise to the non-analytic 
backflow current dominate. Thus, at finite temperature there is a cross-over from the non-analytic result with 
jqp ~ epv s \vfV s /Ao\ for irT <C VfV s -C ttT c to j qp ~ epv s (vfV s /TrT) 2 for < VfV s -C ttT (These expansions are 
discussed in appendix A). Hence, with decreasing surface field, H, the effective penetration depth also crosses over 
from a linear field dependence, SXn ~ A|| (H/H ), for < Ht < H <g; H c , to a quadratic dependence, <5A|| ~ (H/Ht) 2 
for < H < Ht, where the cross-over scale is the field at which the flow energy per excitation becomes comparable to 
the thermal energy, v f v s ~ v f ^H T X ~ nT, or H T ~ (T/A )H . For T/A = 0.001(i.e. T ~ 0.2K and T c = 100 K) 
the cross-over field is Ht — 10 G with H = 10 kG. However, at T = 2 K, the cross-over moves to Ht — 100 G, which 
is a substantial fraction of H c \ ~ 200 G for the cuprates. At temperatures above T ~ 2K the linear field region is 
essentially washed out by the thermal backflow current. Thus, it is essential to work at T <C T c (H c i/Hq) in order 
to minimize the current from the thermally excited quasiparticles. Note that the restrictions on the temperature are 
more severe for observing the linear field dependence, compared to the linear temperature dependence at H = 0, 
because a clean interpretation of the Meissner penetration length requires fields below the vortex nucleation field. If 
vortex nucleation could otherwise be suppressed, then the field range could be extended to H ~ Hq and a linear field 
dependence would be observable over a much larger field range, and correspondingly the restriction on the low-field 
cross-over would be much less severe. In any event we typically assume a field range up to H c \ ~ 250 G in our 
calculations and estimates. 

In Fig. 13 we show the effect of thermal excitations on the velocity dependence of the effective superfluid density 
p s (T, v s ) = js/v s for v s directed along a node. Jl6| ] The intercept shows the usual thermal reduction in p s at v s = 0. 
The cross-over from the linear field dependence at VfV s 3> ttT is clearly seen, and the arrows indicate the value of 
the cross-over velocity, vt = TtT/vf. Figure 14 for the field-dependent Meissner penetration depth contains similar 
information. The T = result is our analytic solution to the nonlinear London equation; for the curves at T ^ we 
converted the current- velocity relations with the same field scaling as we derived for the T — scale, vJvq — > H/ Hq. 
Stojkovic and Vails have recently solved the nonlinear London equations numerically at T ^ 0. p9| Our scaling 
assumption for T ^ agrees well with their numerical solutions for Ht < H c \. Figure 14 gives the estimates of SXn 
for H ~ H c i for several temperatures. The two estimates at T = correspond to H c i/Ho(d.c.) = 0.01, appropriate 
for the d.c. measurement of SXn, and H c \/ Hq(cl.c.) = 0.02, appropriate for the a.c. measurement of 5X\\. Thus, at 
T — > we estimate a change in A|| ~ 30 A in an a.c. measurement for fields directed along the node. At T = OAK 
the cross-over is observable in the calculation, but we expect to see a clear linear field dependence over the field range 
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up to H c \. At T ~ 1 K the resolution is reduced, <5Am ~ 20 A and a linear field dependence is observed over roughly 
60% of the field range. However, by T ~ AK the cross-over field is Ht — H&; the linear field dependence is washed 
out and the change in A|| is less than 15 A. These results are all reduced by a factor of ~ l/y/2 for a field along an 
antinodc. 




Fig. 13 Velocity and temperature dependence of the effective superfluid density, p s = j s /v s , for v s \ \ node in the 
clean limit. The cross-over velocities are indicated by the arrows. 




Fig. 14 Field dependence of the penetration depth for if || node. The cross-over field, Ht = (TVAo)ifo, is 
indicated by an arrow. The verticle marker at 0.01 (0.02) corresponds to H c i/Ho for a d.c. (a.c.) measurement of 
the penetration depth. 



Impurity scattering 

Quasiparticle scattering by impurities also removes the non-analyitc dependence of the current on the condensate 
flow velocity at sufficiently low v s . At T = impurity scattering gives rise to a cross-over velocity v*, or field 
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H* = Hq(v*/vq), that is determined by the energy scale, s* , in the density of states AT(E); above e* <§; VfV s <C 7rT c 
the excitations that are strongly affected by impurity scattering at energies E < e* are only a small fraction of the 
non-thermal backflow current, while at small flow velocities, VfV s < e* the excitations contributing to the non-thermal 
backflow current are strongly modified by impurity scattering. 

The cross-over field scale at T — due to impurity scattering can be obtained from the general expression in cq. 
([52|) for the current; the cross-over velocity is given by VfV* = e* , where e* is the cross-over energy from eq. ([48|). In 
the Born limit (5o <SC 1) the cross-over scale, 

H* = H -i— ~ -H e-^ TA « = -Ho e~^ u ^/^ (Born) , (68) 

is exponentially small for li mp S> £o- Thus, the linear field dependence and anisotropy of the penetration depth will 
be unaffected by impurity scattering in the weak scattering limit. The cross-over field may be much higher, even in 
the dilute limit, for strong scattering. In the unitarity limit (5 — > tt/2) the cross-over field scale is 

H * ~ H ° \ hn (I /r \ -^T (Unitarity) . (69) 
y 21n(A /r„) ^A 

Forr„/A = 10" 4 , which is a good bound for T u obtained from the data of Ref. ( §), we obtain, H* ~ 2.6x10 3 Hq — 
26 G. Thus, even in the unitarity limit there is a large field range, 25 G ~ H* < H < H c \ ~ 250 G in which the 
linear field dependence of <5A|| (H) is expected to hold. Of course it is important to be at low temperatures in order to 
avoid the thermal cross-over. Figure 15 shows our numerical results for the field-dependent penetration depth (d.c.) 
calculated for T/T c — 0.004 and unitarity scattering. The field range is approximately < H < H cl ~ 250 G. In 
the clean limit (r„ = 0) the curvature at very low fields is due to the thermal cross-over discussed earlier. As the 
impurity concentration increases the curvature sets in at higher fields; the arrows indicate the impurity cross-over 
fields calculated from eq. (|69|); note the calculated cross-over field accurately reflects the field dependence obtained 
from the numerical calculations. The results for r u /Ao ~ 10 -4 suggest that a nonlinear Meissner effect with a linear 
field dependence over ~ 80 % of the field range from zero to H cl should be observable in single crystals of comparable 
purity to those of Ref. ( 111) . 




Fig. 15 Field-dependence of the penetration depth with unitarity scattering. The cross-over field H * is indicated 
by an arrow. 



V. ANISOTROPY OF THE IN-PLANE MAGNETIC TORQUE 



Another test of the presence of nodal lines associated with a d x 2 v i order parameter would be to measure the 
magnetic anisotropy energy, or magnetic torque, for in-plane fields. p5| Consider a velocity field v s = v x > x' + ly y' 
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that is not directed along a node or antinode. At T — the projections of the velocity along the nodal lines £' and 
y' give rise to two backflow jets of different magnitudes, 



Jq P = - eN f v j 



Vo J \ V ° J 



(70) 



The important feature is that the current is not parallel to the velocity field, except for the special directions along 
the nodes or antinodes. As a consequence the magnetic field in the screening layer, b, is not parallel to the applied 
surface field, H. This implies that there is an in-plane magnetic torque which acts to align the nodes of the gap, and 
therefore the crystal axes, with the surface field. 

The magnetic free energy of the superconductor, in the presence of the surface field H , is given by 

U = - { d 3 x [ M(H') ■ dH' . (71) 



JVfUm JO 

The integration is carried out assuming that the orientation of H' is fixed; H' = H' (s'm(d)x' — cos(9)y'), and M(H) 
is the equilibrium magnetization for the given value of H , which is easily be found by solving the nonlinear London 
equation in the film geometry (see appendix B). The resulting anisotropic contribution to the magnetic energy is 
obtained by integrating over the volume of the film,^] 

U an {6) = -tLjL AXll ($) [ S in 3 (#) +cos 3 (0)] , (72) 

where A is the surface area of the film, and (see appendix B) 

In the thin and thick film limits we obtain, 

(to) = 1 3 (A) , d< ^\\ (74*) 
[ ' [I >d»A„. 1 ' 

Note that the anisotropy energy is minimized for field directions along the nodal lines, and is maximum for fields 
along the antinodes. 

At finite temperature an analytic expression corresponding to eq. (|70|) is not available except in special limits ( see 
appendix A). In order to calculate the magnetic torque at finite temperature, we note that for |w/Ua| Ao, we can 
write j sx ,y = Psx,yVsx,y (we drop the primes on x and y, but emphasize that the axes refer to two orthogonal nodal 
directions) with 

^ = l-g (VfV s x,y/^A ) , (75) 

where g(x) is a dimensionless function (see the subsections on 'thermal excitations' and 'impurity scattering' in section 
IV). At T — g(x) = \x\, manifesting the non-analytic behavior of the current, while at finite temperature g(x) has 
a linear region for large x (vfV s /T 3> 1), and crosses over to a quadratic region for small x (vfV s /T -C 1). 
The flow velocity and field distribution are determined by the nonlinear London equation, 



-^r 1 - y s x, y {1-5 (vfv s Xiy /fiA a )} = , (76) 



where Am is the penetration length in the limit of zero field. The solution for the velocity field can be obtained by 
performing a perturbation expansion in VfV s /[iAo. For a superconductor occupying the half space z > 0, the solution 



4 The authors of Ref. ( @) erred in calculating the magnetic anisotropy energy by a factor of 3. The corrected result is given 
byeq.@. 
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dz' 



U>x,y 



(77) 



where Q(x) = f x' g(x')dx' . The value of the Ho is determined by the boundary condition, dv s x ,y/dz\ z —o = ±-H y 
thus, 



3 



Hy,? i V 'i H ° ) 



(78) 



\H 



The magnetic torque can be obtained from the magnetic anisotropy energy, t z = —dU an /dO, or equivalcntly M x H. 

The integral of b — V x A can be expressed in terms of v s at the surface, since J °° b x ^ y dz — ±|u s y,x|z=o- Thus, for 
a thick film 









'5)1 








V 3 n 0/ _ 



(79) 



the torque can be calculated from a simple integration of the current-velocity relation. At T = the torque equation 
((79|) reduces to 



1 

t z = — \—H 2 A sin 9 cos 9 (sin 6 — cos 6) , 0<6 l <vr/2, 

47T ilQ 



(80) 



in agreement with the derivative of the anisotropy energy in eq. (|72|) . 

For H = 400 G, A = (2,000/im) 2 and A = 1,400A, the zero temperature maximum magnetic torque r 2 ~ 
(l/12V3ir)H 2 (H/H )A\ - lO" 5 dyne - cm/rad. In Fig. 16 we show results for the magnitude of the torque at 
finite temperature. The torque is only weakly reduced for T < 1 K, but drops rapidly above T ~ 1 K. 




Fig. 16 Magnitude of torque as a function of 9. Note that 9 = 0,n/2 correspond to node positions, 9 — n/4 
corresponds to an antinode position, and the maximum torque (for T = 0) occurs at 9 = ■§ sin" 1 (2/3) ~ 21°. The 
torque has four fold symmetry and points to the nodal positions. 



In order to maximize the torque it is desirable to suppress vortex nucleation so that the torque measurement can 
be performed at a higher magnetic field. In this respect, thin films with dimensions d < A might be desirable, because 
the vortex nucleation field is increased by roughly (X/d) in a thin film. The optimum geometry might be a superlattice 
of superconducting/normal layers with an S-layer thickness, £ <C d < X. In this case the field at each SN interface is 
essentially the external field, and the anisotropy energy is enhanced by the number of S-layer s. 
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Measurements of the transverse magnetization 



In a recent paper Buan, et al. reported |l7j measurements of the transverse magnetic moment induced by an 
in-plane surface field in an untwinncd single-crystal of LuBa2Cu307_,j. The surface field was rotated in the a-b plane 
and the Fourier component of M±(9) cx cos(4#) was extracted and compared with predictions based on the theory of 
Ref. ( Buan, et al. report a measured transverse magnetization signal (~ cos(4#)) of M^ p ~ 0.8 x 10~ 6 emu 

at H = 300 G and T = 2 K, and a resolution limit of 0.3 x 10 -6 emu. Buan, et al. also report a theoretical value 
of Mf eorv ~ 2 x 10- 6 emu, for the same temperature and field, obtained from numerical solution of the London 
equation with the nonlinear current- velocity equation from Ref. ( p!B[). This estimate is 2.5 times the experimental 
signal and nearly an order of magnitude above the resolution limit. The authors conclude that the cos (AO) signal is 
too small to be consistent with a pure d-wave pairing state, that there are no nodes in the gap and that the cos(4#) 
signal is consistent with a higher harmonic of the cos(0) and cos(2(9) signals associated with the shape anisotropy of 
the crystal. 

While it may be that the observed cos(4#) signal is associated with extrinsic effects of geometry, the conclusion that 
the measurement rules out a pure d-wave state with nodes relies principally on Mj_ eorv 3> M^ p . The estimate of 

M theory „ g ^ 1Q _ 6 

emu at T — 2K and H = 300 G is based on the following parameters: A c = 1 \im, A|| = 1, 700 A 
and geometric parameters for the crystal, a = 1.2 mm, b — 0.9 mm and c = 0.07 mm. The transverse magnetic 
moment per unit area is given by M± = t z / (Area H) , which is proportional to H 2 X^ from eq. (80) for T — > 0. Buan, 
et al. argue that finite size effects require that A|| in the formula for M± be replaced by an effective penetration depth 
A ~ 4, 000 A in order to account for c-axis currents. This procedure leads to an increase in the in-plane transverse 
magnetization due to current flow along the c-axis, which is opposite to what is expected. In particular, in a geometry 
with an aspect ratio, c/a ~ 0.07, and a field lying in the ab-plane, say along the 6-axis, the current flows predominantly 
along the a-axis. The 'return current' at the edges flows mainly along the c-axis. The main effect of the c-axis currents 
is a reduction in the area with current flow in the a-b plane, A e ff ~ A(l — 2A c /a). Thus, the in-plane, transverse 
magnetization for the semi-infinite geometry will correspondingly be reduced by roughly A e ff /A = (1 — 2A c /a). The 
reduction factor is tiny for the geometry of Ref. ( jL7|), A e ff/A ~ 0.999; however, the theoretical value for M 1 ^" ^ 
reported in Ref. ( |T^]) is overestimated by a factor of (A/A||) 2 ~ 5.5. Dividing the theoretical estimate of 2x 10~ 6 emu 

by 5.5 gives M^ ^ ~ 0.36 x 10~ 6 emu, very near the resolution limit and below the observed signal at H = 300 G. 

Thus, in our opinion this null result is inconclusive and does not force one to eliminate a pure d-wave state as a 
possible candidate for the order parameter of the cuprates. Experiments designed to minimize shape anisotropics at 
temperatures well below T ~ 2 K should be able to detect the intrinsic anisotropy associated with nodal excitations, 
should they exist. 
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APPENDIX A - TEMPERATURE AND FIELD DEPENDENCE OF THE CURRENT FOR A 

SUPERCONDUCTOR WITH LINE NODES 



In general we need to solve two problems, (i) the self-consistent calculation of the order parameter, A(p/), in 
the presence of a condensate flow field at finite temperature, and (ii) the computation of the temperature and field 
dependence of the current, j s , for a given A(p/;u s ,T). These problems are in general coupled because the flow field 
leads to suppression of the gap parameter, A, which contributes to the field-dependence of the current. We show 
that at low temperatures, T/T c <C 1, and low fields, VfV s /T c -C 1, the contributions to the current from the gap 
suppression can be neglected. 

We also discuss the functional form of the current at low temperatures and low velocity in the clean limit. In the 
limit T/T c <C 1 and VfV s /T <C 1, an expansion in vjv s is valid, and as a result, the leading field dependence of the 
superfluid density, p s , is quadratic in the flow field; while in the limit T/T c <C 1 and VfV s /T 3> 1, the quasiparticle 
contribution to the current is non-analytic. For simplicity we confine ourselves to cases where v s is parallel to a node 
or an antinode, in which case j s and v s are parallel so we drop the vector symbols. 

In the limit |i>/i> s | «T< A we can expand j s in a power series in v s vt, 



j s = -eNfVf v s 



1 



7i 



T 
A^ 



72" 



(vfV a 



TA 



(81) 



where A = A (v s , T) , and_7i , 72 are numerical coefficients of order one whose values depends on the angular slope 
parameter fj, defined in eq. (|45|) . The 71 term is the usual linear dependence of p s on T. The 72 term, being proportional 
to 1/T, signals the breakdown of the Taylor expansion in v s for sufficiently low temperatures. 
The magnitude of the gap has the expansion, 



A (u s ,r) = A (0,0) 



1 — Oil 



T 



A o (0,0) 



o 2 



VfVs 

A o (0,0); VA o (0,0) 



T 



(82) 



Thus, it is clear that the leading correction to j s from the flow field v s is given entirely by the 72 term; gap suppression 
is a higher-order correction. 

In the limit T <C \vfV s \ -C A , we obtain, 



j s = -eNfVf v s 



1 



7i 



, \ V f V s 



A 



72 



T 2 



\vfV s \Ac 



(83) 



The 71 term was calculated earlier, and the 7 2 term can be obtained from the Sommerfeld expansion of integrals 
involving the Fermi function, with \vfV s \ playing the role of the chemical potential for the quasiparticles. Note that 
both coefficents depend on the direction of v s . Also not that the 7 2 term gives a non-analytic current independent of 
\v s \, but the result is only valid for T -C |i>/i> s |. The gap in the above formula is 



A (v Sl T) = A o (0,0) 



l-oi( 



\v f i 



A o (0,0)' 



T 



{VfVs 

A o (0,0)y VAo(0,0) 



(84) 



which gives higher-order corrections to j s . The quadratic correction to j s is unaffected by the gap suppression; the 
lowest-order temperature correction is given by the term cx 7 2 ( \ Vf v s \A ) ■ 



APPENDIX B - NONLINEAR LONDON EQUATIONS 

Conventional superconductor with an isotropic gap 

We include solutions of the nonlinear London equations for a parallel surface magnetic field penetrating into a 
superconductor. Starting from Ampere's equation 

Vx(Vxl) = -J Sl (85) 
c 

consider first an isotropic superconductor with a conventional gap. The constitutive equation for the current is given 
by eq.(|58|). In the absence of phase vortices, V x Vx = 0, we obtain the nonlinear London equation, 
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= 0. 



(86) 



where we have chosen the gauge, V • v s = 0. This gauge choice is consistent with current conservation provided that 
the current is transverse to the spatial variation of the field, which requires that there be no vortices present. Consider 
the geometry in which the superconductor occupies the half-space z > 0. A surface field H = Hx is parallel to the 
interface, and the screening current v s is then parallel to y. The field in the superconductor is given by b = -(dv s /dz)x. 
Continuity of the parallel field at the interface imposes the boundary condition 



dv s 
dz 



y=n — -H . 



(87) 



The field, and therefore the screening current, also vanish deep inside the superconductor. We introduce a dimension- 
less velocity, u = ^/av s /v C7 and distance, C = Z /^(T)', the differential equation for it(£) becomes, 



d 2 u 



with the boundary conditions, 



du , 

d( 



!C=o 



d( 2 

e ^/avfHX _ 



u(l -u 2 ) = 0, 



'o 



u{( oo) = 0. 



(88) 



(89) 



A first integral is obtained by multiplying by du/dC,, integrating, and applying the asymptotic boundary condition, 
u — > as z — * oo, to obtain, 



du 1 

— = — u\ 1 u z . 

d( V 2 

In the physically relevant limit, |u| <§C 1, we replace (1 — ^u 2 ) 1 / 2 - 

Uq 



(1 



iw 2 ), and integrate to obtain, 



V^/4+(l-^/4)e 2 C 



The dimensionless magnetic field is given by 

du 



HQ 



u (l -«o/4)e 



2C 



d( [u 2 /4+(l- U 2 /4)e 2 C] 3/2 



(90) 



(91) 



(92) 



and the constant uq is determined by the boundary condition on magnetic field at the interface, eq. (p9[) ; ho 



-Uq(1 — Uq/4), or for weak nonlinearity, uq 



-ho(l + /Iq/4), which yields the solution, 
h e 2 < 



[h 2 /4+(l-h 2 /4)e 2 i] 



3/2 



(93) 



Note that if we neglect the h\ terms in the denominator we recover the linear response solution to the London equation, 
h = h e - ''. 

A strong surface field produces a reduction in the screening current at the surface and a correspondingly longer 
initial penetration depth. We take the initial decay rate to define the effective penetration depth, 
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(94) 



This definition is equivalent to identifying the effective penetration depth with the surface impedance, i.e. 
l/X e ff(ho) oc j s (0)/H. Note that the effective penetration depth increases with field as expected. In physical 
units 



, , r-Vf e f-H 
h a = Xy/a—-H = \Jct — 
Ac Hq 



(95) 



where Hq = cA/eXvf ~ </>o/A^ — H c and H c is the thermodynamic critial field. 
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d x 2_ y 2 superconductors 

Now consider a superconductor with a d x 2_ y 2 order parameter at T = 0. Choose a coordinate system in which 
the nodes of the gap are directed along the ±x and ±y axes. The nonlinear London equations for the corresponding 
projections of the condensate velocity are 

d 2 u t 

-Ui{l-Ui} = i = x,y, (96) 

where Uj = Vi/vo, C = z /^\\> an d the velocity and length scales are vo = pAg/vf and A = A||. These equations are to 
be solved subject to the boundary condition H = V x A| z= o, which becomes, 

— |c=o = -fto cos 61 , -^|c=o = -«o sin0, (97) 

where 8 is the angle of H measured relative to the node along —y, and ho = |A_ff/u . Note that h ~ H/Hq. The 
differential equation can be solved perturbatively; the first integral is 

du t (. 2 xl/2 

which can be integrated to give 



dc-- Ui [ 1 -r i ) ' (98) 



«i(0 = f -j— t , (99) 

(1 - |u i0 ) cosh(C) + y 1 - |m»o sinh(C)) 



where u^o is the value of the velocity field at £ = 0, which is fixed by the boundary conditions on the field. In the 
limit ho <C 1 we obtain 

u x o = ho cos0 ^1 + —ho cosfc^ , (100) 

u y o = ho sin# ^1 + i/i sin O^j . (101) 

The magnetic field in the screening layer, b = V x A, is given by 

1 f du 



b x /H = - — ^-^j ~+sin0e^ ^1 + -ft o sin0[l - e"^]J , (102) 

b y /H = +±- (jtj^J ~- cose e- c ^l + ^/io cos 0[1 - e" c f) . (103) 

Note that the field in the screening layer is not parallel to the applied surface field. This leads to an in-plane magnetic 
torque acting on the superconductor which tends to align the nodes of the order parameter and the surface field. 

Field penetration in a thin film d x 2_ y 2 superconductor 

A geometry in which the torque anisotropy may be measured is a film of thickness d ~ A||. The field H is oriented 
as in the half-space geometry, parallel to the CuO planes and the surfaces of the film. Choose the origin at the center 
of the film; the boundary conditions for both interfaces are now, 

du du 

-^r\±d/2\ = -h cos9 , _J^| ±d/2A = -h sinO. (104) 
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We solve the differential equation perturbatively. Let Ui(Q — Lj(C) + cti((), where £i(C) is the solution to the 
boundary- value problem for the linearized differential equations, i.e. 

Li(C) =a { sinh(C), (105) 

with a,i fixed by the boundary conditions, 

a x cosh(^-) = — h a cos0 , a y cosh( 7 ^-) = — h a sin . (106) 
2 a 2 A 

The perturbation satisfies the inhomogeneous equation, 

cPoti 2 

-—^--a i = -L i i = x,y. (107) 

The solutions are obtained by writing a{Q — H(()<fr((), where 7i(C) is a solution to the homogeneous equation. The 
function </)(() then satisfies, 

^ + ^^' = -L 2 /H. (108) 



The first-order equation for is obtained by multiplying eq.(108) by the integrating factor H 2 , 

0,(0 = nW Jo L2{x) n{x)dx ' (109) 

The boundary conditions are satisfied at both interfaces by choosing the homogeneous solution to be, TL — sinh(C), 
in which case 

^(0 = ^2 / sinli 3 ^) dx = -y [ cosh (0 - sech 2 (C/2)] . (110) 



One additional integration yields 



0(C) - 4>o = -— [sinh(C) - 2 tanh(C/2)] , (111) 



where cf>o is also chosen such that the full nonlinear solution satisfies the boundary conditions. The resulting solutions 
for the magnetic field are 

b x /H^-^- (^j =+ sine [/5(C) + h sin 6$(Q] , (112) 
by/H = +■! (^) =- cos 0[/3(C) + ^ocos 0$(0] , (113) 



/3(C)=cosh(C)/cosh(^), (114) 



$(C) = 3sech d (— ) x | C osh(C)sinh^( — )cosh( — ) - cosh( — ) sinh J (|) cosh(|) } . (I IT,) 

Note that the nonlinear correction is largest at a distance of order A/2 from the interface. 
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